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SUMMARY 

The  Dynamic  Engine  ANalysIs  program,  DEAN,  Is  a  FORTRAN  code  Implemented 
on  the  IBM/370  mainframe  at  NASA  Lewis  Research  Center  for  digital  simulation 
of  turbofan  engine  dynamics.  DEAN  Is  an  Interactive  program  which  allows  the 
user  to  simulate  engine  subsystems  as  well  as  a  full  engine  systems  with  rela¬ 
tive  ease.  The  non-linear  first  order  ordinary  differential  equations  which 
define  the  engine  model  may  be  solved  by  one  of  four  Integration  schemes,  a 
second  order  Runge-Rutta,  a  fourth  order  Runge-Kutta,  an  Adams  Predictor- 
Corrector,  or  Gear's  method  for  stiff  systems.  The  numerical  data  generated 
by  the  model  equations  are  displayed  at  specified  Intervals  between  which  the 
user  may  choose  to  modify  various  parameters  affecting  the  model  equations  and 
transient  execution.  Following  the  transient  run,  versatile  graphics  capa¬ 
bilities  allow  close  examination  of  the  data. 

DEAN'S  modeling  procedure  and  capabilities  are  demonstrated  In  this  paper 
by  generating  a  model  of  a  simple  compressor  rig. 


INTRODUCTION 

System  simulation  Is  a  powerful  tool  In  the  aerospace  Industry.  These 
simulations  can  be  us -J  In  conjunction  with  testing  of  a  specific  system  In  an 
attempt  to  anticipate  possible  problems  and  reduce  risk.  Also  simulations  can 
be  utilized  In  subsystem  development,  specifically  In  the  design  and  develop¬ 
ment  of  controls  for  turbofan  engines. 

These  simulations  must  have  several  attributes  In  order  to  be  useful 
tools.  They  must  be  reasonably  accurate  and  computationally  efficient. 

Currently  there  are  several  programs  In  use  that  simulate  turbofan 
engines.  NNEP  (ref.  1)  Is  a  generalized  code  for  steady-state  analysis  of 
arbitrary  engine  configurations.  Another  code,  DYNGEN  (ref.  2)  simulates  the 
dynamics  of  specific  engine  configurations.  A  complex  hybrid  (analog/digital) 
computer  simulation,  HYDES  (ref.  3)  simulates  dual  shaft  turbofan  engines  with 
up  to  three  air  streams.  These  simulations  are  limited  as  to  their  complexity 


♦Currently  with  NASA  Lewis  Research  Center. 


because  the  hybrid  computer  has  hardware  and  memory  limitations.  DIGTEM 
(ref.  4),  a  recently  written  code  allows  steady-state  as  well  as  transient 
calculations  to  be  performed.  OIGTEM  has  the  flexibility  to  model  various 
engine  configurations  which  are  subsets  of  a  two-spool,  two  stream  turbofan. 
DIGTEM's  Interactive  capability  Is  limited  and  turbofan  engine  components  can¬ 
not  be  developed  separately.  None  of  the  above  packages  have  the  extensive 
flexible  graphics  provided  In  DEAN. 

The  Dynamic  Engine  ANalysIs  program,  DEAN  has  been  written  specifically 
to  model  two-spool,  two-stream  turbofan  engines  and  turbofan  engine  components. 
However,  DEAN  can  be  easily  reconfigured  to  model  other  dynamic  systems  such 
as  turboshaft  engines,  Internal  combustion  engines,  or  electrical  power 
systems.  User  Interaction,  modular  engine  component  routines,  various  Inte¬ 
gration  methods,  and  extensive  graphics  have  been  Incorporated  Into  the  soft¬ 
ware  to  enhance  the  simulation  capability  of  DEAN. 

Implemented  under  the  TSS  operating  system  on  an  IBM/370  mainframe  com¬ 
puter,  DEAN  was  written  with  the  Intent  of  simulating  dynamic  systems  Inter¬ 
actively.  The  program  utilizes  the  IBM/370's  Interactive  debugging 
capabilities,  graphics  libraries,  and  the  IMSL  FORTRAN  library  which  allows 
the  user  to  develop  models  on-line  In  the  sense  of  finding  modeling  errors  and 
Initiating  transients.  The  Interactive  capability  also  allows  the  user  to 
analyze  system  response  to  external  Inputs  such  as  perturbations  In  fuel  flow, 
nozzle  area,  and  bleed  flows.  Temperatures,  pressures,  flows,  and  other  model¬ 
ing  parameters  otherwise  known  as  state  variables  are  displayed  continuously 
during  the  simulation.  Also,  part  of  DEAN'S  user  Interaction  Includes  Its  user 
friendliness  which  Is  achieved  by  using  menus  with  logical  default  selections 
during  program  execution. 

There  Is  a  natural  division  of  components  In  a  turbofan  engine,  l.e., 
Inlet,  fan,  compressor,  combustor,  high  pressure  turbine,  low  pressure  turbine, 
afterburner,  nozzle,  bypass  duct,  and  engine  control.  DEAN  Is  designed  to 
simulate  a  turbofan  engine  with  a  set  of  modular  routines,  one  for  each  engine 
component.  The  user  may  use  the  provided  component  models  for  his  simulation 
or  he  may  choose  to  provide/develop  his  own  model  for  any  or  all  of  the 
components . 

Generally  the  models  are  based  on  first  principles,  l.e.,  conservation  of 
mass,  momentum,  and  energy,  and  can  be  formulated  as  a  system  of  first  order 
nonlinear  differential  equations.  DEAN  uses  these  equations  and  one  of  several 
numerical  Integration  techniques  to  approximate  the  system  response.  The 
Integration  technique  may  be  selected  by  the  user  which  best  approximates  the 
system  being  simulated.  The  appropriate  use  of  the  available  Integration 
schemes  will  be  discussed  later. 

DEAN  supports  an  extensive  graphics  capability  for  displaying  results  of 
a  simulation.  This  capability  Includes  the  ability  to  plot  any  parameter. 
Including  time,  versus  any  other  parameter,  superposition  of  performance  maps 
on  a  plot,  calculation  of  new  parameters  based  on  existing  parameters,  the 
option  of  having  one  or  two  plots  per  page,  extensive  labeling  Including  date, 
time,  and  confidentiality  of  data,  ability  to  plot  the  data  for  time  sub¬ 
intervals,  and  the  capability  to  save  the  data  for  later  analysis.  The  plots 
generated  by  DEAN  are  of  report  quality. 


There  are  several  support  packages  required  to  model  turbofan  engines. 
Generally  the  engine  performance  maps  are  bivariate  In  nature,  l.e.,  compressor 
pressure  ratio  Is  a  function  of  corrected  speed  and  corrected  flow.  Also  var¬ 
ious  thermodynamic  properties  are  required  such  as  Cp  (air  and  fuel),  Cv, 
and  y,  the  ratio  of  specific  heats,  all  of  which  are  single  valued  functions. 
Therefore  subroutines  have  been  supplied  to  Interpolate  single  and  bivariate 
maps.  Also  there  are  subroutines  for  reading  free  formated  datasets  and  a  sub¬ 
routine  Is  Included  for  simple  proportlonal/lntegral  control. 


DEAN  STRUCTURE 

DEAN  was  designed  to  be  a  modular  program.  Subroutines  were  written  to 
do  specific  tasks  such  as  prompting  the  user  for  a  steady-state  point  or  run¬ 
ning  a  transient  simulation.  This  manner  yields  several  advantages.  DEAN  can 
be  expanded  and  modified  easily  because  only  specific  modules  may  need  minor 
changes.  Also,  the  modularity  In  the  program  expedited  the  software  develop¬ 
ment  phase.  Further,  the  programming  flow  of  DEAN  Is  more  easily  understood. 
Thus  DEAN  Is  a  highly  structured  program  which  allows  the  user  to  easily  Iden¬ 
tify  and  modify  components. 

DEAN  consists  of  46  subroutines  that  are  called  by  one  main  program. 

Since  there  are  so  many  subroutines,  only  the  more  Important  ones  will  be  dis¬ 
cussed  here.  Figure  1  shows  a  simplified  block  diagram  of  the  subroutine 
Interaction  In  DEAN. 


Figure  1.  -  Simplified  diagram  of  Dean  subroutines. 
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The  main  program  (DEAN)  Is  responsible  for  controlling  the  entire  program. 
It  presents  the  user  with  the  following  options:  (1)  system  definition,  (2) 
map  definition,  (3)  steady-state  operating  point  calculations,  (4)  transient 
simulation,  (5)  data  analysis  (graphics),  (6)  termination  of  DEAN  -  and  prompts 
for  a  response.  When  the  user  has  entered  an  appropriate  response,  the  program 
calls  the  subroutine  associated  with  that  response. 

Subroutine  DEFSYS  Is  prompts  the  user  for  the  components  of  the  system  to 
be  simulated.  By  using  a  menu  of  components  (table  1.),  DEFSYS  allows  the  user 
to  define  a  system  for  simulation.  The  order  of  these  components  corresponds 
to  the  configuration  of  a  turbofan  engine.  If  the  user  desires  to  simulate  a 
system  other  than  a  turbofan  engine,  he  must  redesignate  the  components  to 
correspond  to  his  system.  For  example,  when  simulating  a  turboshaft  engine, 
the  low  pressure  turbine  component  would  contain  a  model  for  torque  loading. 
With  minimal  software  modifications  to  various  menus  In  DEAN,  the  system  con¬ 
figuration  may  be  formally  adapted.  The  user  will  also  be  prompted  for  the 
required  model  conditions  such  as  Inlet  temperature,  pressure,  and  fuel  flow 
If  the  combustor  and/or  the  afterburner  are  being  simulated  without  a  fuel 
control . 

From  the  components  chosen  by  the  user,  subroutine  MAPSET  decides  which 
bivariate  maps  may  be  needed  for  the  simulation  and  prompts  the  user  to  execute 
one  of  two  tasks  related  to  Initialization  of  these  maps.  The  rough  map  data, 
supplied  by  the  user,  resides  In  specified  format  In  a  dataset.  When  beginning 
a  simulation  MAPSET  processes  the  rough  data  Into  a  uniform  grid  for  the  bivar¬ 
iate  Interpolation  routines  and,  upon  request,  rewrites  the  data  over  the 
original  dataset.  The  processing  of  rough  data  needs  to  be  done  only  once 
because  MAPSET  can  determine  when  the  data  In  the  dataset  Is  rough  data  or  the 
pre-processed  grid.  DEAN  Is  capable  of  processing  eight  bivariate  maps  which 
have  been  designated  as  shown  In  table  2.  The  user  may  wish  to  change  the 
designated  maps  and,  may  do  so  as  long  as  no  more  than  eight  bivariate  maps 
are  specified.  Also,  upon  request,  MAPSET  will  graphically  display  the 
processed  map  data  as  each  map  Is  processed.  This  allows  the  user  to  verify 
the  maps  specified  by  the  data  he  has  supplied. 


TABLE  1.  -  LIST  OF  COMPONENTS 
SIMULATED  BY  DEAN 


Component  number 

Turbofan  engine  component 

1 

Inlet 

2 

Fan 

3 

Compressor 

4 

Combustor 

5 

High  pressure  turbine 

6 

Low  pressure  turbine 

7 

Afterburner 

8 

Nozzle 

9 

Bypass 

10 

Rotor  dynamics 

11 

Main  fuel  control 
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TABLE  2.  -  AVAILABLE  BIVARIATE  HAPS 
SUPPORTED  BY  DEAN 


Hap  number 


Description 

Fan  pressure  ratio 
Fan  temperature  ratio 
Compressor  pressure  ratio 
Compressor  temperature  ratio 
High  pressure  turbine  pressure  ratio 
High  pressure  turbine  enthalpy 
Low  pressure  turbine  pressure  ratio 
Low  pressure  turbine  enthalpy 


Subroutine  SSOPPT  Is  responsible  for  providing  the  simulation  with  an 
Initial  starting  point  or  steady-state  point.  There  are  two  mechanisms  the 
user  may  choose  to  determine  such  a  point.  Generally  this  point  consists  of  a 
vector  of  state  variables  that  represent  the  pressure,  temperatures,  flows, 
shaft  speeds,  etc.  of  the  system.  There  Is  one  state  variable  for  each  differ¬ 
ential  equation  In  the  system  model  as  well  as  any  parameters  calculated  by 
empirical  or  theoretical  equations  such  as  choked  flow  for  example.  The  first 
mechanism  Is  to  Input  the  states  directly  from  the  keyboard.  DEAN  will  prompt 
the  user  for  the  states,  component  by  component  as  specified  by  the  system 
definition  portion  of  DEAN.  An  option  Is  provided  to  use  this  Input  as  an  Ini¬ 
tial  guess  In  conjunction  with  the  governing  equations  and  Iterate  to  a  steady- 
state  point  by  a  modified  Newton- Raphson  technique.  An  alternate  method  for 
determining  a  steady  point  Is  to  use  the  Initial  guess  as  a  starting  point  for 
a  transient  and  run  the  transient  until  the  dynamics  settle  out.  Once  a  steady 
state  point  has  been  calculated.  It  may  be  stored  with  the  system  configuration 
for  future  reference.  This  allows  the  user  to  select  a  saved  steady-state 
point  rather  than  entering  the  states  each  time  a  simulation  Is  run. 

The  subroutine  TRNSNT  Is  responsible  for  running  the  transient  simulation. 
TRNSNT  prompts  the  user  for  one  of  four  Integration  schemes  and,  after  an 
acceptable  Input,  requests  more  specific  Information  such  as  Initial  time  step, 
length  of  transient  time,  frequency  at  which  states  should  be  displayed  at  the 
terminal,  and  also  whether  the  user  wishes  to  Interact  directly  with  the  sim¬ 
ulation.  Oata  generated  by  the  transient  simulation  Is  stored  at  specified 
Intervals  for  graphical  analysis  at  a  later  time. 

The  subroutine  PLOT  Is  responsible  for  generating  graphics  output  from 
OEAN.  This  subroutine  supports  various  options  highlighted  by  table  3.  These 
options  are  discussed  briefly  later  In  the  paper. 

The  last  group  of  Important  subroutines  are  listed  In  table  4.  These  sub¬ 
routines  contain  all  of  the  modeling  equations  required  by  the  simulation. 
Included  In  these  subroutines  are  data  for  the  component  geometry,  data  for 
single  variable  functions,  code  to  store  data  for  graphical  analysis,  bivariate 
map  calls,  evaluation  of  specific  thermodynamic  properties  such  as 
Cn  =  f (temperature) ,  etc.  These  subroutines  can  be  written  by  the  user  who 
wishes  to  develop  an  analytical  model,  as  opposed  to  using  those  provided. 
Currently  OEAN  Is  set  up  to  Integrate  up  to  one  hundred  state  variables. 


TABLE  3.  -  GRAPHICS  PACKAGE  OPTIONS 
FROM  SUBROUTINE  PLOT 


Option 

Description 

1 

Read  plot  labels  from  labels  dataset 

2 

Read  plot  data  from  dataset 

3 

Save  plot  data  on  dataset 

4 

Define  or  modify  plot  parameters 

5 

Display  a  menu  of  available  plot  parameters 

6 

Change  the  time  Interval  of  the  data  plotted 

7 

Automatically  make  a  hardcopy  of  each  plot 

8 

Generate  a  standard  cartesian  plot 

9 

Superimpose  a  map  on  a  standard  cartesian  plot 

10 

Put  data  classification  labels  on  all  plots 

TABLE  4.  -  LIST  OF  DEAN'S 
COMPONENT  SUBROUTINES 


Subroutine 

Description 

EQUNS 

Subroutine  for  "global"  equations 

INLET 

Inlet  duct  model 

FAN 

Fan  model 

COMP 

Compressor  model 

BURNER 

Combustor  model 

HPTURB 

High  pressure  turbine  model 

LPTURB 

Low  pressure  turbine  model 

AFTBRN 

Afterburner  model 

NOZZLE 

Nozzle  model 

BYPASS 

Bypass  duct  model 

MFC 

Engine  control  model 

ROTDYN 

Dynamic  rotor  balance  equations 

It  may  be  noted,  that  while  the  component  subroutine  names  correspond  to 
those  of  a  turbofan  engine,  they  may  also  be  used  to  model  components  of  vir¬ 
tually  any  dynamic  system. 


INTERACTIVE  ASPECTS  OF  DEAN 

The  simulation  program,  DEAN,  has  five  major  areas  of  user  Interaction. 

The  first  area  Involves  defining  the  system  components  (listed  In  table  1). 

OEAN  will  prompt  the  user  for  the  types  of  components  required  for  the  Simula 
tlon,  and  thus  subsystems  of  the  turbofan  engine  may  be  modeled,  l.e.,  compres¬ 
sor  and  nozzle.  Based  on  the  component  models  used,  OEAN  will  prompt  the  user 
for  specific  Information  such  as  Inlet  pressure  and  temperature  of  the  model 
or  fuel  flow  for  a  combustor. 

The  second  area  of  Interaction  Is  map  grid  definition  for  bivariate  maps 
(z  =  f[x,y]).  DEAN  will  prompt  the  user  for  Information  concerning  the  defini¬ 
tion  of  the  grid.  The  program  can  presently  handle  as  many  as  eight  bivariate 


6 


maps  as  listed  In  table  2.  Such  maps  are  usually  required  for  fan,  compressor, 
high  pressure  turbine,  and/or  low  pressure  turbine  performance  data. 


The  next  area  of  Interaction  concerns  the  Initial  values  of  the  states  for 
the  simulation.  The  user  must  Initially  enter  values  for  the  states  which  give 
the  Integration  schemes  a  starting  point.  These  values,  which  usually  corre¬ 
spond  to  a  steady  point  are  entered  from  the  keyboard.  The  user  may  also 
choose  to  select  a  previously  entered  steady  point. 

The  most  useful  area  of  Interaction  Is  during  the  actual  system  simula¬ 
tion.  As  the  transient  proceeds,  the  program  will  save  data  on  specified 
parameters  (up  to  100)  with  a  maximum  of  1000  points  per  parameter.  The  method 
of  saving  the  data  Is  an  Integral  part  of  the  component  subroutines.  During 
the  transient  the  user  may  change  various  parameters  based  on  the  simulation's 
output. 

The  final  area  of  Interaction  Is  In  data  analysis.  DEAN  supports  a  menu 
driven  graphics  package  (table  4)  with  capabilities  for  basic  mathematical 
manipulation  of  the  transient  data. 


MODULAR  ENGINE  COMPONENTS 

A  turbofan  engine  Is  generally  comprised  of  those  components  listed  In 
table  1.  In  the  Dynamic  Engine  ANalysIs  program  there  Is  a  subroutine  allo¬ 
cated  for  each  of  these  components.  The  subroutines  consist  of  FORTRAN  code 
for  evaluating  derivatives  and  storing  data  for  analysis.  Also,  there  Is  code 
to  generate  parameters,  perform  map  lookups  for  single  and  bivariate  maps,  and 
evaluate  thermodynamic  properties  for  specific  components.  These  component 
subroutines  are  always  called  by  DEAN  In  the  order  shown  In  table  1. 

This  modular  simulation  structure  allows  the  user  to  develop  turbofan 
engine  systems  and  subsystems,  component  by  component.  This  method  Is  useful 
because  It  allows  easier  debugging  of  both  coding  and  theoretical  problems. 
Currently,  the  only  restriction  Is  that  any  subsystem  must  have  Its  components 
ordered  as  per  table  1.  For  example,  a  turbofan  core  consisting  of  a  compres¬ 
sor,  a  combustor,  and  a  turbine  In  line  may  be  simulated.  However,  a  change 
In  component  arrangement  or  system  topology  would  require  minor  code  changes 
and  possible  renaming  or  addition  of  component  subroutines. 


INTEGRATION  SCHEMES 

The  mathematical  model  of  the  system  Is  represented  generally  by  ordinary 
first  order  nonlinear  differential  equations.  There  are  several  methods  which 
can  be  used  to  Integrate  these  equations.  DEAN  supports  the  following  Inte¬ 
gration  schemes:  (1)  second  order  Runge-Kutta,  (2)  fourth  order  Runge-Kutta, 
(3)  Adams  predictor-corrector,  and  (4)  Gear's  method  for  stiff  systems.  Each 
of  these  Integration  schemes  has  Its  own  appropriate  applications. 

The  second  and  fourth  order  Runge-Kutta  schemes  (ref.  5)  execute  simi¬ 
larly.  The  second  order  scheme  requires  fewer  function  evaluations,  and  thus, 
less  CPU  time.  However,  the  fourth  order  scheme  Is  more  accurate.  In  general, 
these  explicit  methods  should  be  used  when  model  frequencies  and  accuracy 
requirements  are  high. 
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The  Adams  predictor-corrector  Integration  (ref.  6)  has  several  appealing 
attributes.  It  Is  a  variable  order  (up  to  12),  variable  time  step  method.  The 
algorithm  Is  a  multlstep  process  In  which  there  Is  a  predictor  for  extrapolat¬ 
ing  to  the  solution  and  a  corrector  for  Interpolating  a  solution  based  on  the 
extrapolation.  The  Interpolation  procedure  Involves  Iteration  and,  therefore, 
requires  more  computation  time  than  the  one-step  Runge-Kutta  (fourth  order 
method)  If  there  are  more  than  two  Iterations  In  the  corrector  process. 

The  Gear's  algorithm  for  Integrating  differential  equations  (ref.  7)  was 
developed  for  systems  that  have  a  wide  range  of  frequencies  (stiff  systems). 

A  turbofan  engine  can  be  considered  as  such  a  system.  The  Gear's  method  Is  a 
backward  difference  scheme  capable  of  variable  order  (up  to  five)  and  variable 
time  step.  It  Is  useful  for  simulations  that  require  a  substantial  amount  of 
computer  time  for  function  evaluations. 


GRAPHICS  CAPABILITIES 

The  last  of  DEAN's  main  menu  options  Is  the  Data  Analysis  package.  This 
set  of  subroutines  allows  the  user  to  manipulate,  as  well  as  plot,  data  gen¬ 
erated  during  a  transient.  The  main  subroutine  which  controls  execution  of 
this  set  of  code  Is  subroutine  PLOT.  PLOT  contains  a  menu  of  ten  options  for 
the  user  to  choose  from. 

The  ten  PLOT  options  allow  the  user  to:  (1)  read  user  specified  plot 
labels  corresponding  to  the  parameters  saved  during  a  transient,  (2)  read  plot 
data  stored  In  a  dataset  from  a  previous  transient,  (3)  write  transient  data 
Into  a  dataset  for  future  reference,  (4)  perform  basic  mathematical  manipula¬ 
tion  of  the  transient  data  and/or  modify  plot  title  and  labels,  (5)  display  a 
list  of  the  plot  labels,  (6)  change  time  Interval  over  which  the  data  Is  plot¬ 
ted,  (7)  automatically  hardcopy  plots  (assuming  hardcopy  capability  exists), 

(8)  produce  cartesian  plots  of  specified  parameters,  (9)  plot  specified  para¬ 
meters  with  a  bivariate  map  (corresponding  to  table  2)  super-imposed  on  the 
data,  (10)  write  security  classification  labels  on  plots  of  the  transient  data. 


CONCLUDING  REMARKS 

DEAN  Is  a  versatile  simulation  program  written  In  FORTRAN  IV.  Designed 
for  Interactive  dynamic  simulation  of  turbofan  engine  systems  and  subsystems, 
DEAN  has  successfully  simulated  a  number  of  component  models,  such  as  the  com¬ 
pressor  rig  example  In  the  Appendix  and  a  full  turbofan  engine.  Further,  the 
Interactive  graphics  capabilities  have  proven  to  be  useful  for  simulation 
development  and  analysis. 

DEAN  Is  supported  by  Interactive  debugging,  the  IMSL  library  (9th 
edition),  and  a  FORTRAN  accessible  graphics  package. 


APPENDIX  A 


COMPRESSOR  RIG  MODEL  EXAMPLE 

The  Dynamic  Engine  Malyhs  program  can  best  be  understood  by  actually 
developing  a  simulation  of  a  system.  Here,  a  simple  simulation  of  a  compressor 
rig  Is  used  to  demonstrate  the  capabilities  of  DEAN.  Also,  the  model  develop¬ 
ment  will  not  overshadow  the  use  of  OEAN  as  a  simulation  program. 

First  the  model  Is  defined.  Let  us  assume  that  the  compressor  rig  can  be 
modeled  as  a  constant  speed,  one-dimensional,  multi-lumped  volume,  Incompress¬ 
ible  system  with  three  major  components:  (1)  Inlet  volume,  (2)  compressor,  and 
(3)  turbine  nozzle.  Note  that  these  three  elements  correspond  to  those  In 
table  1.  The  Inlet  will  be  modeled  as  a  single  lumped  volume  and  the  compres¬ 
sor  will  also  be  represented  as  a  single  volume  employing  the  actuator  disk 
concept  suggested  by  Greltzer  (ref.  8)  for  compressor  stall.  The  turbine 
nozzle  will  be  modeled  by  using  a  choked  flow  relation  between  temperature  and 
pressure. 


Generally  three  parameters  provide  communication  between  components:  here 
(l)pressure  (P),  (2)  mass  flow  ($),  and  (3)  temperature  (T)  are  used. 

These  parameters  are  calculated  by  Integrating  the  non-linear  first  order  dif¬ 
ferential  equations  that  result  from  the  conservation  of  mass,  momentum,  and 
energy  and  the  equation  of  state  In  differential  form.  Figure  2  shows  a  typi¬ 
cal  lumped  volume  with  the  various  parameter  designations.  The  general  form 
of  these  equations  for  such  a  lumped  volume  analysis  Is  as  follows,  where  p 
Is  density,  t  Is  time,  V0]  Is  volume,  R  Is  gas  constant  of  air,  L  Is  the 
volume  length,  A  Is  cross  sectional  area,  Fext  denotes  sum  of  frictional 
and  pumping  forces,  Q  denotes  the  net  heat  transfer  to  the  volume,  W  denotes 
net  work  done  by  the  volume,  Cv  Is  the  constant  volume  specific  heat  of  the 
fluid,  and  y  the  ratio  of  specific  heats. 
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Figure  2.  -  Typical  lumped  volume  diagram. 


Figure  3.  -  System  block  diagram. 


TABLE  5.  -  STEPS  REQUIREO  TO  DEVELOP  A  SIMULATION 


Step  number 

Instructions 

1 

2 

3 

4 

5 

6 

Define  the  system  to  be  modeled 

Olvlde  the  system  Into  components 

Write  dynamic  equations  for  each  component 
Write  FORTRAN  component  subroutines 

Create  datasets  to  be  used  by  simulation 
Access  and  execute  DEAN 

From  these  equations,  models  of  an  Inlet  and  compressor  can  be  developed. 

Table  5  shows  the  general  steps  Involved  In  developing  the  simulation.  A  block 
diagram  of  the  system  Is  shown  In  figure  3.  Using  the  block  diagram  and  the 
general  conservation  equations  the  Inlet  equations  may  be  written  as  follows: 

dP,  P_T,  P,  dT_ 

pressure  (state)  gT  =  7—  -  *c>  *  j~  dT 


dW  gc  (  .2 

mass  flow  (momentum)  dT  =  L~  )P1A1  ~  P2A2  “  R  W2 


dT2  RT? 

temperature  (energy)  =  7-7 -  "  V  "  *<TyT2  "  V 

2  0 1 , 2  ( 

where  Fext  Is  modeled  as  Fext  =  -RlSl2  (R  =  constant  to  give  appro¬ 
priate  pressure  drop).  Also,  for  the  Inlet,  P  and  T  are  constant  Inlet 
conditions.  The  compressor  Is  modeled  In  much  the  same  way.  The  equations  for 
the  compressor  are  as  follows: 

dP_  RT.  P,  dT 

pressure  (state)  =  7-^-  (*c  -  *n>  -  fT  dT 


dWc  9C 

mass  flow  (momentum)  ^7-  =  —  jp2A2  “  P3A3  +  Pext 
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where  the  values  of  P2>map  and  P3(map  are  the  pressures  corresponding 
to  a  steady  state  compressor  performance  map  In  terms  of  pressure  ratio  (pres¬ 
sure  ratio  Is  a  function  of  corrected  flow  and  corrected  speed)  and  T2  map 
and  T3 jinap  are  temperatures  corresponding  to  another  steady-state  perform¬ 
ance  map  In  terms  of  temperature  ratio  (temperature  ratio  Is  a  function  of 
corrected  flow  and  corrected  speed.)  The  turbine  nozzle  will  be  modeled  by 
using  the  following  choked  flow  nozzle  equation: 


mass  flow 


W  =  a  P  I-*-1 — 
n  rr3  h  ♦  1 


RT3  1/2 


y  ♦  1 

Y  ♦  1 


The  governing  equations  are  now  coded  Into  separate  subroutines  for  each 
component.  Once  these  subroutines  are  written  for  all  of  the  component  models 
the  user  must  collect  all  the  data  required  for  the  model.  For  the  compressor 
example  there  are  two  data  requirements;  (1)  compressor  map  data  and  (2)  geom¬ 
etry  data. 


Once  all  of  the  preliminary  work  Is  done  for  the  simulation  DEAN  may  be 
run.  A  typical  run  Is  similar  to  that  discussed  In  the  INTERACTIVE  ASPECTS  OF 
DEAN  section. 


For  this  particular  example,  the  nozzle  exit  area  Is  ramped  closed  during 
the  transient  while  compressor  shaft  speed  Is  held  constant  as  shown  In 
figure  4.  The  Adam's  method  was  chosen  as  the  Integration  scheme  with  an 
Initial  time  step  of  0.0001.  The  length  of  the  transient  was  0.50  sec  which 
required  6.63  sec  of  CPU  to  Integrate  six  states.  This  disturbance  causes  a 
compressor  Instability  called  surge.  Figure  5  shows  graphically  the  Instabil¬ 
ity  superimposed  on  the  appropriate  performance  map.  Figures  6  through  8  show 
time  traces  of  Inlet  and  compressor  exit  flow,  temperature,  and  pressure, 
respectively. 
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Figure  6.  -  Time  traces  of  weight  flow. 
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Figure  7.  -  Time  traces  of  weight  flow  and  temperature. 
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